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\—\ ! Abstract 

We use a modified version of the Nambu-Jona-Lasinio model which implements 
the QCD trace anomaly to calculate the two photons decay width of the glueball 
(/o(1500)) and of the two scalar mesons (/o(1370), /j=o(1710)) to which it is mixed. 
We investigate the effect of this mixing over the coupling constants of the /o states 
to the quarks. 
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1 Introduction 



The scalar resonance /o(1500) which has been discovered at LEAR by the 
Crystal Barrel Collaboration in the reaction pp (at rest) — > 3ir° , nrj-n- [1] and 
nrj V [2], was immediately considered as a good candidate for the scalar 
glueball. Indeed, the observed mass ((1520±25) MeV in [1] and (1545±25) 
MeV in [2]) was near the value predicted by the lattice-QCD calculations [3] 
of that time. It is also in agreement with the analysis of recent data on the 
decay /o(1500) — > KK [4] which however requires a mixture between the 
glueball and the neighbouring qq states [5]. Another signal that strengthens 
this identification lies in the fact that, in the radiative decay J/ip — > 7 + 4-7T, 
a sharp peak is seen in the 4n spectrum, having quantum numbers ++ and 
mass and width also compatible with / (1500) [6]. 
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The mixing between the glueball and the qq states has also been noted by 
Sexton et al. [7]. They however identify the nearly glueball state with /j(1710) 
whose mixing with / (1370) and / (1500) requires that J = [5], [8]. 

In Refs. [8], [9], it has been claimed that good tests of the glue and qq contents 
of the glueball would consist in studying its partial strong decay widths into 
7i7i , K K , r)T) ; , r)f] . We plan to calculate them in a forthcoming paper [10]. 

In the present paper, we follow the idea of Close et al. [5] who estimate the 
two photons decay widths of the glueball and of the mesons to which it is 
mixed. Emission of photons is surely a good test to investigate the glueball or 
gg-nature of the /o states since gluons do not couple to photons. The authors 
suggest a mixing between / (1370), / (1500) and / (1710) for which they 
study two different schemes, according to the preceding remarks. Firstly, they 
consider the case where /o(1500) exhibits large contents in glue and in qq (u 
and s) excitations while / (1710) has its largest component in the ss chanel. 
This allows to reproduce the /o(1500) — > KK data [4]. The second scheme 
consists in considering the glueball lying above the ss member of the nonet 
[7]: /o(1710) can then be the glueball while /o(1500) is mainly a ss excitation. 
Within these two schemes, Close et al. [5] have estimated the relative strength 
of the 27 widths for the three /o states: 

/o(1370) : / (1500) : / (1710) « 12 : 1 : 3, (scheme 1) 

/o(1370) : / (1500) : / (1710) « 13 : 0.2 : 3. (scheme 2) 

Whatever the mixing, the decay width of /o(1500) is always the smallest. 
Experimental estimation of the 27 widths would then be a good test of the 
general idea of qq and glue mixing. Indeed, if the width of the / (1500) was 
found to be larger than one of the others, the idea of mixing would be de- 
stroyed. Note that the decay / (1370) — > 27 is the only one that has been 
measured experimentally : r_f ( 1370 ) = 5.4 ± 2.3 KeV [11]. 

The model we use to calculate the two photons decay widths is a modified 
version of the SU(3) NJL model [12] which implements the trace anomaly of 
QCD [13]. It is called the scaled NJL model. It entails the introduction of a 
scalar field \ whose mean value xo can be identified with the vacuum gluon 
condensate [14]. This x field couples to the quarks so that our model can 
describe the processes schematically shown in Fig. 1. In our model, it is the 
vacuum gluon condensate that fixes the mixing [15]: the smaller xo, the larger 
the mixing. 

However, in the domain of Xo which allows to roughly reproduce the meson 
masses, / (1710) remains mainly a ss excitation. Its coupling to the up quark 
is weak. In contrast, the strength of the coupling to the u or s quark is roughly 
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Fig. 1. Schematic 27 decay of the glueball and of a scalar state. 

of the same order of magnitude for / (1370) and for / (1500). In our model, the 
three states have glue content. / (1500) is the one which yields to pure glueball 
if no mixing and it is the state that always has the largest glue content even 
if that of /o(1370) can become important. We then work in a scheme similar 
to the scheme 1 presented above. 

The parameters of the model are fixed to reproduce the pion and kaon masses, 
the pion decay constant and the masses of ao(1450) [ud,ud, -jz(uu — dd)\ and 
/o(1500). The mass of the largest member of the nonet (/o(1710)) will then 
vary as well as that of /o(1370). The latter is very sensitive to the value of the 
vacuum gluon condensate while the former keeps a value not far away from 
1710 MeV, whatever xo ■ 

For completeness, let us add that the 27 decay of the scalar mesons has also 
been studied in Ref. [16] within a relativistic quark model with linear con- 
finement and instanton-induced interaction. / (1500) and /o(980) are there 
considered as mixing between uu and ss states. Another relativistic treatment 
using a OGE model has been presented in Ref. [17]. 

Our paper is organized as follows. Section 2 recalls the useful tools of the scaled 
NJL model. Section 3 is devoted to the derivation of the decay widths of the 
mesons and of the glueball. We insist on the fact that the mixing between 
the three fields generating the glueball and the scalar mesons modifies the 
transition amplitudes. For instance, /o(1370) can decay into 27 through the 
production of a ss pair. The estimation of the widths assumess the knowledge 
of the masses of the mesons and the glueball as well as their coupling constants 
to the quarks. These quantities are defined in Section 4. We discuss our results 
in Section 5. Finally, Section 6 gives our conclusions. 



2 The model 



The model used in this paper is described in Refs. [12,18] under the name of 
" A-scaling model" . We recall here some of the useful tools for the understand- 
ing of the present work. We start from the vacuum SU(3) effective Lagrangian: 



£eff(q, q, x) = q(-id^u + m )q 

which yields to the bosonized effective Euclidean action 
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+ / <i 4 x— — V9 a v? a + / d^xC x (2) 



once one has integrated out the quark degrees of freedom. The meson fields 
write: 

V?a = (<r a ,7T a ), r a = (A a ,i7 5 A a ), a = 0,...,8 (3) 



where the A a are the usual Gell-Mann matrices with Ao = a/2/31. We choose 
to work in the isospin symmetry limit and the quantity m stands for the 
diagonal matrix diag(m u , m u ,m s ). The trace anomaly of QCD is modelized 
using a scalar dilaton field x that is intimately related to the gluon condensate 

X oc (Gl) 1/2 [12]: 
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£x = o(<W + T^Vln-V - (X 4 - X G ))- (4) 



Since we are only interested in the scalar sector, the axial anomaly which 
would give the 77-?/ mass difference is not considered here. Our model contains 
six parameters: the current quark masses (m u ,m s ), the strengths (a 2 ,b 2 ), the 
gluon parameter xg an d the cut-off A introduced to regularize the diverging 
quark loop. Four of these parameters are adjusted to reproduce the pion mass 
(m,,-), the weak pion decay constant (/V), the kaon mass (m^) and that of 
/o(1500). We then have two free parameters that we choose to be the vacuum 
gluon condensate Xo (related to xg [18]) an d the constituent up quark mass 
M u . The latter is connected to the strength a 2 using the usual gap equation 
[12] which minimizes the action. Note that the strange constituent quark mass 
M s and xo obey similar stationary conditions. M u will have a large value (725 
MeV) in order to reproduce ao(1450) [8] and xo wn l remain free. 



3 The decay widths 



The widths for the two photons decay of the scalar isoscalar mesons are given 
by 
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for i = / (1370), /o(1710) and /o(1500). In Eq. (5), to, represents the mass of 
the meson, Q its quadrimomentum while k\ and k 2 stand for the quadrimo- 
menta of the emerging photons. The mesons are assumed to be at rest. In a 
pure NLJ, the transition amplitude for /o(1500), identified with the glueball, is 
identically zero; rj ( 1370 ) only involves the propagator of the up quark (in the 
limit of isospin symmetry) while T f (yj W \ is expressed in terms of the strange 
quark propagator. 
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Fig. 2. Feynman direct diagrams for 27 decay of /o(1370) (a) and of /o(1710) (b). 
The direct diagram (Fig. 2. a) writes 

T /o(137 o) = N c (e 2 u + e 2 d ) gf {im E^k^MU (6) 

with 



u = J -nSi^D {{i 2 + M u )^ 3 + M u ) lv {ix + M u )} 
1 

x 



{ql + MlM + MiM + MD 

while the direct diagram (Fig. 2.b) has the similar form: 

2/o(ino) = N c e 2 s g% {17w ft(hK(k 2 )S 



(7) 



with 



S = I ~^%^D {(i 2 + M s ) 7m (4 + M s ) lv (i x + M s )} 

J ( 27 

1 

X fe 2 + M2)(g| + M2)(gl + M2)- 



(9) 



In the scaled NLJ, the mesons /o(1370) and /o(1710) are not anymore pure 
uu or pure ss excitations. Equations (6) and (8) have then to be modified in 
order to take these new effects into account: 



r /o(13 ro) =N c (e 2 u + f&gffignriihyMU 

+ N c elg}*l m e;(h)E:(k 2 )S 



(10) 



T M mo) = N c eig} s s {1710) e;(k 1 )et(k2)S 



+ N c (e 2 u + e 2 d )g*^ 1710) e;(fa)e:(fa)U. 
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The quark meson coupling constants g* uu , g* 



/o(1370), /o(1710)) differ 



from the ones introduced in Eqs. (6) and (8). Their analytical expressions will 
be established in Section. 4. 

Having some content in qq excitations the glueball /o(1500) can now decay: 



T mm =N c (e 2 u + e^g^flihK^U 

+ N c e*g}°l m e;(h)e:(k 2 )S. 

Summing over the final polarization of the photon and using Q 
qi = q - fa, <?2 = q + fa, q 3 = q, one gets : 
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with 
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and a = e 2 /47r the fine structure constant^. We will work with a large con- 
stituent up quark mass, M u = 725 MeV, in order to reproduce do (1450). 
However, even with such a large mass, the mesons /o(1710) and the glueball 
(1500) always lie above the unphysical pair uu creation threshold 2M U . In the 
same way, /o(1710), which aquires a mass larger than 2M S = 1728 MeV, is 
always above the ss threshold. The function J u should then be calculated in 
the complex plan: 

J u ^J u (Q>-iO+)\ Q2= )2 , (17) 



where Tj is the width of the meson for its uu decay and with a similar ex- 
pression for J s . This width as well as the mass m^ can be calculated from a 
Dyson- Schwinger equation (see below) which amounts to two coupled equa- 
tions for Tj and m^. The latter cannot be decoupled without a suitable and 
model-dependent description. In the same way, expression (17) is here approx- 
imated by | J u (— rrii 2 — i0 + )|. Note that the results are not significantly affected 
by the prescription used above the threshold. However, it would worth using 
a confining model to get rid of the unphysical decay channel. Indeed, Celenza 
et al. [19] have shown that a modification of the scalar vertex by a confining 
nonlocal interaction increases the mass of the a of ~ 250 MeV. This feature 
could surely modify quantitatively the results we present here. 



4 Meson masses and coupling constants 



Equations (13-15) involve the masses of the mesons as well as their coupling 
constants to the quarks u and s. We now show how to calculate them, insisting 
on the way the coupling constants are modified due to the coupling between 
the various fields. The effective action (2) can be expanded up to second order 
in the fluctuating parts of the meson fields (do, erg) and of the dilaton field x '■ 

1 „ (°o,-9 

I {2) {a ,a 8 ,x) = ^H(^o, g , cr 8 , ff , x,)S _1 v 8 ,- q 
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3 Detailed calculations for ttq — ► 77 can be found in Ref. [18] 



S x is a 3x3 matrix with 9 nonvanishing elements that can be read in Ref. 
[12]. It can be diagonalized, yielding to 
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with 



diag A u (q 2 ) =V 1 (q 2 )S- 1 V(q 2 ), 



(20) 



and V the q 2 dependent eigenvector matrix. Expression (19) can be approxi- 
mated by expanding the elements of the diagonal matrix around their respec- 
tive zero, that is around the respective meson masses: 
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In Eq. (21), we have used the fact that 

diag Ku(-m 2 ) 



im-q- 



(21) 



(22) 



and we have defined [15] 

G- 2 = 9 g2 diag \Ki{-m 2 )] = \v T (-m 2 )(d q2 S- l )V(-m 2 ) 



(23) 



where the eigenvector matrix V has its first row calculated at — m 2 

-m 2 f /-.Tini and the last one at —ml = 



-m 



its second row at —m\ 
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The new physical fields i associated with the three /o states write: 



G7% = G; 1 J2v ij 1 



u ' m 1) a i 



/o(1500)- 



(24) 



with o"i = cr , o- 2 = o- 8 , 0-3 = %. 

The interaction term in the quark loop (In) between the quarks and the me- 
son fields (see Eq.(2)) reflects an interaction term in the Lagrangian of the 



form q(<Jo\o + agAs)?. Using Eq. (24), the latter becomes an interaction term 
between the quarks and the <f> i : 

q(a \o + o"8As)g = uAiu + dA 2 d + sA 3 s (25) 



with (there is no summation over the index i) 

Ai= E ^((^) lt V lj + (\s) lt V 2j )G J i = l,2,3 (26) 

j=l,2,3 



showing that the two mesons as well as the glueball are indeed coupled to the 
three quarks with the respective coupling constants: 

9*r = 9f=[(*o)nVij + (As)nVy Gj (27) 

^=^0)33^ + ^8)33^]^, (28) 

j = l(/o(1370)),2(/ (1710)),3(/o(1500)). 

In a pure SU(3) NJL, without mixing with the dilaton, Eq. (26) reduces to 

A i = $' i G i i = l,2 (29) 

so that Eq. (23) defines the coupling constant of the mesons to the quark and 
4>\ ((fil) 1S onr y coupled to the quark u(s). 

In the present approach, we will use Eqs. (22), (27), (28) to calculate the 
mass of the mesons and their respective coupling to the quarks. Note that it is 
the nonvanishing value of <7f"nsoo) anc ^ #/o(i50o) ^ na ^ a ^ ows the glueball decay 
/ (1500)^77- 

As mentioned above, some caution must be taken when the uu or ss thresholds 
are reached. Due to the mixing with the glueball, it is not possible to solve 
analytically diag [Ajj(g 2 )] (see Eq. (22)). In order to illustrate the threshold 
problem, we will drop here this mixing. The diagonalization then yields: 

A* (g 2 ) = AN c F u { q 2 W + 4M 2 ) + a 2 X 2 ^ (30) 

771 

A 22 (g 2 ) = 4iV c F s (g 2 )(g 2 + 4M 2 ) + a\ 2 ^- (31) 

with 

FM 2 )=ld A k T -^ — — i=- — — T i = u,s. (32) 



(k-ff + M? (k+ff + M? 



Whenever the mesons lie below the threshold their mass is given by: 



A 11 (— ^/ O (i37o)) — 0, A 22 (— TOj ( 1710 )) — 0. 



(33) 



When the threshold is reached, Eq. (33) has to be modified to take into account 
the fact that the meson acquires a width due to its nonphysical decay into uu 
or ss. This width and its mass are solutions of : 
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or more explicitly: 

AN C F U \-(m fo +ie- iT fo ) 2 } -(m fo - iT /o ) 2 + 4M 2 






2 2 '""U n /oc\ 

a XoTT- = °> ( 35 ) 



with similar expression for / (1710) with u replaced by s. In order to get 
rid of the unphysical Tf , we perform here calculation with the approximated 
equations [15] 



4V, 



-(m fo +iey\ (-m} o +4M u J ) + a 2 Xo^ = 0, i = u,s. (36) 



Other approximated expressions for the right-hand side of Eq. (35) can be 
found in the literature [20]. To our knowledge, Eq.(35) has never been solved 
exactly. 



5 Results 



The Xo behavior for the mass of the three fo states is given in Fig. 3. For 
large values of \oi the states take the mass they would have without coupling 



1^/0(1370) 



2M,. 



1450 MeV 



m /o(mo) 



2M, 



1728 MeV in the chiral 



limit). The mass of /o(1710) is quite stable while an acceptable value for the 
mass of /o(1370) limits the value of xo to a domain between 300 and 450 MeV 
(see Table 1). 

The same kind of behavior appears at the level of the coupling constants of the 
f states to the quarks u (Fig. 4, Eq. (27)) and s (Fig. 5, Eq. (28)). Note that 
Fig. 4 and Fig. 5 exhibit the modulus of the coupling constants. Whenever 
the phase factor (mj o — 4Mf) is negative, the corresponding coupling constant 
also becomes negative in such a way that the two contributions u and s add up 
(see Eqs. (13)— (15)). Here again, g* uu and g* ss are stable for the / (1710) with 
g* ss an order of magnitude larger than g* uu indicating that /o(1710) is nearly 
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Fig. 3. xo behavior of the mass of the three states /o(1370), /o(1500) and /o(1710). 

a pure ss state. For large Xo, ffffino) - *^ in contrast, the coupling constants 
g*uu anc j g*ss are Q £ ^e same order of magnitude for / (1370) and for / (1500). 
At large xo, <?/ofi37o) S oes t° zero, reflecting the fact /o(1370) is a pure uu state 
when the mixing is turned off. Finally, /o(1500)— »27 is forbidden in our model 
for large xo, illustrated by the asymptotic behavior of ^"nsoo) anc ^ fl'/ofisoo) 
that vanish when xo~ ¥C, °- However, the decrease of the coupling constants of 

Table 1 

Masses and widths in 2 7 of the hybrids / (1370), /o(1710) and / (1500). 







masses (MeV) 


r(eV) 


r / r /o(1710) 






/o(1370) 


1341 


348 


51.9 


Xo 


= 400 MeV 


/o(1500) 


1500 


102 


15.2 






/o(1710) 


1753 


6.7 


1 






/o(1370) 


1305 


617 


69.7 


Xo 


= 350 MeV 


/o(1500) 


1500 


112 


12.7 






/o(1710) 


1755 


8.85 


1 






/o(1370) 


1250 


1258 


102.3 


Xo 


= 300 MeV 


/o(1500) 


1500 


123 


10.0 






/o(1710) 


1757 


12.3 


1 
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/o(1500) is slower than the others. 
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Fig. 4. xo behavior of the coupling constant to the quark u, g* uu , of the three states 
/o(1370), /o(1500) and / (1710). 
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Fig. 5. xo behavior of the coupling constant to the quark s, g* ss , of the three states 
/o(1370), / (1500) and / (1710). 

The masses and the coupling constants so obtained allow to estimate the 
27 decay widths (Eqs. (5),(13)-(15)). Results are shown in Fig. 6. The three 
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widths decrease with increasing xo since for xo ~^ oo, at least one of the corre- 
sponding coupling constants vanishes. Fj ( 1710 ) is always the smallest because 
of two additionnal effects: in the RHS of Eq. (14) the first term is small due 
to the small value of g*f U n 7 i ) while the second one nearly vanishes due to the 
fact that ra/ (mo) ~ 2M S . For the two other states, the relative amplitude of 
the decays depends on the value of xo- 
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Fig. 6. xo behavior of the width into 27 of the three states /o(1370), /o(1500) and 
/o(1710). 



For small Xo, r /o(i370) > r /o (i 50 o) 
Cisco) and mJ o(1370) - AMf 



due to four cumulative effects: 
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m- t f, „™ - 4M 2 



2 

/o(1500) 

Decomes smaller whi 



with % = u,s. For large xo, 
etheRHSofEq. (15) keeps 



two nonvanishing terms up to quite large values of xo- One then has r^ ( 1370 ) < 
r/ (i5oo)- Note that in the chiral limit, the three widths go asymptotically to 
zero. (The factors (m 2 — 4M 2 ), which lead to a vanishing of the decay widths 
of /o(1370) and / (1710) in the pure NJL model (no mixing) in the chiral 
limit, are a consequence of the chiral symmetry in the scalar sector. This was 
already shown by Bajc et al. Ref. [21].) 



One then sees that our model provides meson and glueball widths that can vary 
from one to three orders of magnitude. Moreover, their relative magnitude is 



/„ (1370) /T/ (1710) 



also xo dependent (see Fig. 7). Taking r/ ( m o) as reference, T 

varies of two orders of magnitude and rj ( 1500 )/rj ( m0 ), while more stable, 

can however vary of a factor 5 in the considered domain of xo- 



We have indicated in Table 1 the value of the widths as well as their ratio to 
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for three typical values of xo compatible with the masses of the /o 
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Fig. 7. xo behavior of the ratios r /o(1370) /r /o(mo) and r /o(1500) /T /o(mo) . 
states. 

Our results are in complete desagreement with the statement of Close et al. [5] 
presented in the introduction: the relative strengths are not reproduced, even 
qualitatively since the width r_f ( 150 o) is always larger than rj ( 1710 ). Moreover, 
the width of / (1370) is always too small. To get the lower bound of the 
experimental value (5.4 ± 2.3 keV), one needs to consider gluon condensates 
as small as 270 MeV yielding a mass of ~ 1200 MeV and a still higher ratio 
r/r /o (i7i ). 

In order to compare our results with those of Ref. [16] it seems appropriate 
to compare our ratio r /o (i 37 o)/T /o (i5oo) with their result for r /o ( 98 o)/r /o (i5 9 o). 
Indeed, in both cases the coupling of the two mesons to the uu and ss channels 
is important. Of course our model provides an additionnal glue content. One 
finds a quantity which varies from 3 to 10 while their corresponding result lies 
between 6 to 11. 



6 Summary and conclusions 



The model we developed in Refs. [14,15] implements the QCD trace anomaly 
of QCD by the introduction of a dilaton scalar field \- (Other models of 
this type are on the market. See for example [22,23] and references therein.) 
Since the NJL model is not renormalizable, a cut-off must be introduced to 
regularize the diverging integrals. This cut-off breaks the scale invariance of 
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the quark loop that can be restored by the replacement A — » A%. This entails 
a mixing between the three scalar isoscalar fields. Put in an other way, one 
emerges with three scalar "hybrids" /o which are mixing of glue and uu and 
ss excitations. We followed the idea of Close et al. [5] who identify these 
hybrids with /o(1370), /o(1500) and /o(1710). In our model, one has one free 
parameter, the vacuum gluon condensate xo that fixes the strength of the 
mixing. In the domain of xo (300 MeV < xo < 450 MeV) which reproduces in 
an acceptable way the masses of the states, the mixing is quite large, reflected 
by the large values of the coupling constants to both types of quarks u and s, 
especially for / (1370) and / (1500). The latter is identified with the glueball 
in the sense that it is the state that would yield pure glue if there was no 
mixing. In that case, it could not decay into 27. Here this decay is allowed 
as well as that of the two other states. The relative magnitude of the widths 
largely depends on the value of the gluon condensate. However, whatever xo? 
the width of /o(1710) is always the smallest. Our results are then at variance 
with that of Close et al. [5] according to which it is the width of /o(1500) that 
is always the smallest. They also claimed, that the idea of mixing should be 
revised if it was found experimentally that rj ( 1500 ) could exceed the width 
of the other states. We have however presented a model which yields totally 
different results while still including mixing and have shown that the reduction 
factors (m 2 — AM 2 ) (playing almost no role for the /o(1500)) are the key to 
understand the discrepancy between our results and those of [5]. Of course, 
our model contains some drawbacks, the most important being surely the 
lack of confinement. However the statement of Ref. [5] should be considered 
with some reserve and we suggest that people using model with dilatons [23] 
coupled to scalar states investigate the problem. 
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